Effect of pressure on the Raman modes of antimony 
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The effect of pressure on the zone-center optical phonon modes of antimony in the A7 structure 
has been investigated by Raman spectroscopy. The Ag and Eg frequencies exhibit a pronounced 
softening with increasing pressure, the effect being related to a gradual suppression of the Peierls-like 
distortion of the A7 phase relative to a cubic primitive lattice. Also, both Raman modes broaden 
significantly under pressure. Spectra taken at low temperature indicate that the broadening is 
at least partly caused by phonon-phonon interactions. We also report results of ab initio frozen- 
phonon calculations of the Ag and Eg mode frequencies. Presence of strong anharmonicity is clearly 
apparent in calculated total energy versus atom displacement relations. Pronounced nonlinearities 
in the force versus displacement relations are observed. Structural instabilities of the Sb-A7 phase 
are briefly addressed in the Appendix. 

PACS numbers: PACS: 78.30.-j 63.20.-e 71.15.Nc 62.50.+p 



I. INTRODUCTION 

Like other semi-metals of the group V-A, antimony 
crystallizes at ambient conditions in the trigonal A7 
structure (space group Rim, No. 166). ^ Antimony is 
known to transform at w8.5 GPa to a high-pressure 
phase^ that is referred to as Sb-II in the recent hterature. 
At 28 GPa Sb converts to the body-centered cubic mod- 
ification Sb-IIIi^ Diffraction studies^ show Sb-II to 
adopt a tetragonal self-hosting structure with an incom- 
mensurate guest-host arrangement. In addition, incom- 
mensurate modulations occur in Sb-IIi^ The structural 
sequence for Sb under pressure is further complicated by 
the observation of another incommensurate phase, called 
Sb-IV, in a narrow pressure range from 8 to 9 GPa (up- 
stroke), i.e. intermediate between Sb-I and Sb-II. ^'^ 

Here we are interested mainly in the effect of pres- 
sure on optical phonon modes of the A7 phase of Sb. 
The A7 structure (Fig. [ij is a rhombohedrally distorted 
variant of a cubic primitive structure. The formation of 
paired layers results in a doubling of the number of atoms 
in the primitive rhombohedral cell. In this context, the 
term Peierls distortion is appropriate because the pairing 
is driven by a gain in electronic energy (formation of a 
pseudo-gap in the electronic density of states at the Fermi 
level E-p) which outweighs the cost in Ewald energy >S In 
the following we will mostly refer to the hexagonal de- 
scription of the structure. The non-primitive hexagonal 
unit cell contains six atoms. For conversion relations be- 
tween rhombohedral and hexagonal cells we may refer to 
Ref.|l3or|lll 

Group theory for the A7 structure predicts three zone- 
center optical modes, i.e. a singly degenerate Aig mode 
(totally symmetric) and a doubly degenerate Eg mode. 
The atoms are displaced along the C3 axis in the Aig 
mode and perpendicular to that axis in the Eg mode. 
The Aig mode is the 'Peierls distortion mode'. In the un- 
folded Brillouin zone of the cubic primitive structure, the 



Aig mode corresponds to a longitudinal acoustic mode 
at the R-point, whereas the Eg mode corresponds to a 
transverse acoustic mode at R. 

Previous Raman studies of antimony at ambient pres- 
sure have addressed the resonance effects^i^ the second- 
order scattering by optical modes the effect of tem- 
perature on the first-order modes and the coupling of 
the first-order Raman modes to electronic excitations 
The coherent Aig phonon generation in Sb by short laser 
pulses is the subject of several more recent reports.-- Ra- 
man experiments on Sb under pressure were performed 




FIG. 1: (Color online) The A7 structure of Sb. Two distor- 
tion modes lead from the parent cubic-primitive (cPl) to the 
A7 (hR2) structure: a rhombohedral elongation along a cubic 
[111] direction and a pairing of the original cubic (111) lattice 
planes. Atoms in red mark opposite corners of the distorted 
cube. The resulting layered structure has a 3-1-3 atomic coor- 
dination. The non-primitive hexagonal unit cell is indicated 
in the figure. 
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by Richter et ali^ who measured frequency shifts up to 
0.6 GPa and observed a softenmg of the zone-center op- 
tical modes, not only for Sb, but also the A7 phases of 
As and Bi. Optical phonons in As were studied up to 
30 GPa using Raman spectroscopy-'^^ whereas in the case 
of Bi femtosecond pump and probe techniques were ap- 
plied to explore the effect of pressure up to 3 GPa on 
the modulation of the optical response by coherent Aig 
phonons liSiSS 

In this work, we have investigated the full stability 
range of the A7 phase of Sb using Raman spectroscopy. 
We also report Raman spectra of Sb-II measured up 
to 16 GPa. Our discussion focuses on results for the 
A7 phase, specifically the pressure-induced softening of 
mode frequencies up to the first phase transition and the 
observed strong pressure dependence of the Raman line 
widths. Experimental results for phonon frequencies of 
the A7 phase are compared to frozen-phonon calculations 
performed within the density functional framework using 
optimized structural parameters. A side aspect of our 
calculations concerning structural instabilities of the A7 
phase in the vicinity of the Sb-I to Sb-II transition is 
presented in the Appendix. 



II. RAMAN SPECTROSCOPY 

The Raman scattering measurements were performed 
on crystalline antimony of 99.999% purity. A diamond 
anvil cell (DAG) was used for pressure generation. Ex- 
periments were initially performed with nitrogen as a 
pressure medium and then, in order to avoid the Ra- 
man scattering of solid nitrogen,^ continued with the 
standard 4:1 methanol-ethanol mixture (M/E) as well as 
with helium; the latter ensures the best hydrostatic con- 
ditions in combination with low background. In all ex- 
periments, pressures were measured by the ruby lumines- 
cence method i 'Si-?^ Raman spectra as a function of pres- 
sure were recorded in backscattering geometry (normal of 
the sample surface parallel to the C3 axis) using a triple- 
grating spectrometer (Jobin-Yvon T64000) in combina- 
tion with a liquid-nitrogen-cooled charge-coupled device 
(GGD) detector. The resolution was set to < 2 cm~^. An 
argon ion laser was used for excitation at a wavelength 
of 514.5 nm. The laser was focused to a spot of 50 /xm 
in diameter. The power incident at the sample surface 
was 15-20 mW. At that power level we did not notice any 
indication for laser-induced heating of the sample when 
mounted inside the DAG. 

Figure 121 shows a subset of Raman spectra of antimony 
measured at room temperature and at different pressures 
up to 16.5 GPa. Spectra at 0, 2, and 5.6 GPa are for 
the A7 phase, those at 8.3 and 9.1 GPa are for phase 
mixtures, and those above 10 GPa are for Sb-II. The 
frequencies of the two observed Raman modes of the A7 
phase decrease with increasing pressure. The Aig feature 
is slightly asymmetric at low pressure, and with increas- 
ing pressure a shoulder develops at its high energy side. 



A broad feature between 235 and 320 cm~^ seen at am- 
bient pressure can be attributed to second-order scatter- 
ing by optical phonon modesii^ Under applied pressure, 
that feature could not be separated well enough from the 
background scattering to allow for useful information to 
be obtained. 

Schwarz et al.^ pointed out that the transition to Sb-II 
is sluggish and that a pure phase Sb-II can be obtained 
only above 13(1) GPa. On the other hand, Degtyareva 
et al. obtained a pure phase Sb-II at 9.1(1) GPa (up- 
stroke) . Our Raman spectra from different runs indicate 
a transition regime extending from about 8 to 10 GPa. 

After passing the transition regime, the Raman spectra 
look quite different from those of the A7 phase. At least 
three new Raman features appear as indicated by the ar- 
rows in Fig. 121 AH of them shift to higher frequency with 
increasing pressure. We could not reproducibly identify 
other distinct Raman features of Sb-II, except for a broad 
shoulder on the high energy side of the most intense Ra- 
man peak. 

It has been reported that, at ambient pressure, the 
Eg and Aig Raman features exhibits a Fano-like line 
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FIG. 2: Raman spectra of antimony at different pressures 
(T — 300 K) . Spectra are scaled to about the same maximum 
intensity; the scaling factors are given in the figure. The 
broad feature seen at about 270 cm~^ (see blown-up section 
of the bottom spectrum) is attributed to second-order Ra- 
man scattering of the A7 phase. The spectra at 8.3 GPa and 
9.1 GPa correspond to a phase mixture of A7 and high pres- 
sure phase(s), those at 11.7 and 16.5 GPa to Sb-II. Arrows 
mark Raman features characteristic of Sb-II. 
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FIG. 3: (Color online) Raman peak frequencies of Sb as a 
function of pressure at room temperature. Closed symbols are 
for increasing pressure, open symbols for decreasing pressure. 
Data were collected in runs using different pressure media. 
The transition regime between the A7 phase and phase II is 
indicated by vertical dashed lines. Lines passing through the 
data points represent results of least squares fits (for parame- 
ters see Table EJ. The inset shows the splitting of the optical 
mode frequencies of the A7 phase as a function of pressure. 
There, the line is a guide to the eye. 



shape which is believed to arise from the interaction 
between electronic intraband scattering and one-phonon 
scattering. These asymmetries are weak, however. So, 
for the purpose of extracting peak positions and approx- 
imate line width parameters from the Raman data for 
the A7 phase, a fitting by symmetric pseudo-Voigt line 
shapes was found to be adequate. For Sb-II, peak fre- 
quencies were determined by simply locating intensity 
maxima. The obtained frequencies of Raman features 
of Sb for pressures up to 16 GPa are shown in Fig. O 
Within the scatter of the frequency data, the pressure 
effects are fully reversible. 

The observed ambient-pressure Aig and Eg mode fre- 
quencies are 149.8 cm~^ and 110 cm~^ at room temper- 
ature [values measured for one of the samples at 2 K are 
155.7(7) and 117.0(3) cm^^]. For a pressure of 7.5 GPa, 
the total relative shift is -18% and -24% for the Aig and 
Eg modes. The magnitude of these shifts is similar to 
what has been found for the optical modes of arsenic in 
the more extended stability range (0 - 25 GPa) of its A7 
phase.ifi> The nonlinear uj{P) behavior of the Eg mode of 
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TABLE L Frequencies of the Raman features of the Sb-A7 
and Sb-II phases. The pressure dependence is described by 
uj{P) — LOo+a P+b P^+c where frequency lo is in cm~^ and 
pressure P in GPa. For the Sb-I phase, the zero-pressure fre- 
quencies oJo, the linear pressure coefficients A, the quadratic 
one (B) for the Eg mode, and the cubic coefficient C for the 
Aig mode are given. The zero-pressure mode Griineisen pa- 
rameters 7o are determined from the linear coefficients a and 
the reported bulk modulus Bo = 39.7 GPa.* For the Sb-II 
phase, we list Raman peak frequencies lj at 9.1 GPa and the 
average linear coefficients a for the pressure range from 9 to 
16 GPa. The average values 7 of the mode Griineisen pa- 
rameter given for that range are estimated assuming a bulk 
modulus value of 100 GPa. 
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Sb can be described by the usual combination of a linear 
and quadratic pressure coefficient, that of the Aig mode 
is approximated well by a combination of a linear and a 
cubic term in pressure. Sets of polynomial parameters 
characterizing the pressure dependence of mode frequen- 
cies and mode Griineisen parameters are given in Tabled 

The frequencies of the Raman bands attributed to the 
phase Sb-II are 121, 126, and 135.6 cm~^ at 9.1 GPa. The 
frequencies of all three features exhibit an essentially lin- 
ear pressure dependence up to 16 GPa. The correspond- 
ing mode Griineisen parameters 7 given in Table ^ were 
estimated using a bulk modulus value of Bq = 100 GPa 
as extracted from pressure- volume data given in Ref. 0- 

The structure of Sb-II consists of an incommensu- 
rate guest- host arrangement. For the host structure 
with tetragonal symmetry and space group /4/mcm|3, 
a group theory analysis gives five zone-center Raman- 
active modes of symmetry Aig, Big, B2g, Eg(l), and 
Eg (2). The corresponding displacement patterns are 
those of the aluminum sublattice in CuAl2P^ Only three 
Raman features could be clearly identified in our spectra 
for Sb-II. Not having investigated polarization effects, we 
do not propose an assignment at this point. So, concern- 
ing Sb-II, the main result of the present study is that 
the frequencies of Raman peaks presumably due to opti- 
cal phonons are similar to that of the Aig mode of Sb-I 
and, different from Sb-I, the optical phonons of Sb-II 
harden with increasing pressure and the values of the 
mode Griineisen parameters are larger than one. 

In view of the observed changes in the phonon frequen- 
cies under pressure, we briefly mention two related prop- 
erties of Sb. First, there is the negative slope of the melt- 
ing line of Sb-I. The melting temperature drops from 904 
K at ambient to about 830 K at 5 to 5.7 GPa (the triple 
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FIG. 4: (Color online) Raman line widths (FWHM) of A7 
antimony as a function of pressure at room temperature. Dif- 
ferent shapes of symbols are used for different pressure media. 
Closed symbols are for increasing pressure, open symbols for 
decreasing pressure. Dashed lines are guides to the eye. Note 
that the pressure-induced broadening is fully reversible. Ar- 
rows indicate the change in width upon cooling to 5 K (cf. 
spectra shown in Fig.|KJ. 



point between liquid, Sb-I and high pressure phase) 
Within a 'single-phase' approach to melting and in the 
spirit of the Lindemann melting criterium, a microscopic 
explanation of the negative melting slope would have to 
take into account that at least part of the phonon spec- 
trum softens considerably with increasing pressure. The 
second remark concerns the pressure-induced supercon- 
ductivity of antimony. The superconducting transition 
temperature Tc of Sb-I increases with increasing pres- 
sure, reaching 0.7 K at 8.5 GPai^l At the phase transition 
to Sb-II, Tc jumps to 3.7 K and then it drops continu- 
ously to 2 K at 25 GPa.^^ Qualitatively, the increase of 
Tc in the A7 phase can be understood to arise from an 
increase of the density of states NjEp) (i.e. the filling of 
the pseudogap, cf. Refs. l28l and l2^ and, in view of the 
strong phonon softening, an enhanced electron-phonon 
couplingi^ In the phase Sb-II, N{Ef) is much larger com- 
pared to that of Sb-I at ambient pressure i^SiSS The effect 
of pressure on N{Ep) of Sb-II appears to be not very 
pronouncediSS So, one can speculate that a hardening of 
relevant phonon modes, as indicated by our Raman data 
for optical modes, accounts for the decrease of Tc with 
increasing pressure in Sb-II. 

Besides a large phonon softening in Sb-I, the second 



very pronounced effect of pressure is the change in line 
widths of the Aig and Eg Raman peaks. The pressure 
dependence of the peak widths (FWHM) is shown in 
Fig. 0] Error bars take into account some ambiguity in 
separating a peak from the background and/or shoul- 
ders. Between ambient pressure and 8 GPa, the width 
of the Aig feature increases by about a factor of three, 
from ~6 cm"-'^ to ~17 cm~^. The Eg-related peak starts 
with a width of ~9 cm~^ at ambient and the width in- 
creases substantially to ^^30 cm~^ at 8 GPa. For the 
N2 medium (solid at 2.5 GPa and changing phase at 6 
GPa), the width becomes even larger under pressure; so, 
the observed line width of the Eg feature seems to de- 
pend sensitively on non-hydrostaticity or shear strength 
of the pressure medium. The data obtained with alcohol 
and helium medium are considered to reflect the effect 
of truly hydrostatic stress on the Raman line widths. It 
is important to note that the pressure-induced broaden- 
ings were found to be fully reversible within experimental 
scatter. This also applies to runs which, during upstroke, 
passed into the Sb-II phase regime. 

We have performed one high-pressure experiment at 
low temperature. The sample was loaded in He medium, 
taken to 7.1 GPa at 300 K, and then cooled to 5 K. Dur- 
ing coohng the pressure dropped somewhat to 5.8 GPa 
due to experimental factors. The high-pressure spectra 
at 300 K and 5 K are shown in Fig.|Sl together with room- 
and low-temperature spectra measured at ambient pres- 
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FIG. 5: (Color online) A room-temperature Raman spectrum 
of Sb at 7.1 GPa compared to the spectrum measured after 
cooling to 5 K. Upon cooling, the pressure dropped to 5.8 
GPa. Spectra were measured for a helium pressure medium. 
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sure. At zero pressure and when cooling, the width of the 
Aig drops from 5.4 to 4.1 cm~^, that of the Eg mode from 
8.9 to 6.6 cm^^. The low-temperature widths are larger 
compared to those reported in Ref./** presumably due to 
crystal imperfections induced by the sample preparation 
procedures. At high pressure the change in width is 14 
^8.5 cm-i for Aig and 28 20 cm"! for Eg. The 
latter changes are indicated by arrows in Fig. ^ These 
changes in peak widths demonstrate that about half of 
the phonon decay rate at room temperature is thermally 
activated. Spectral weight components like the shoulder 
on the high-frequency side of Aig also depend on tem- 
perature. 

A possible explanation of the line widths changes un- 
der pressure would be the coupling to acoustic phonons 
via cubic anharmonicity.^^''^*' At ambient pressure, the 
Eg mode frequency represents the lower frequency limit of 
the optical phonon branches, separated by about 10 cm~^ 
from the acoustical mode regimejiS Based on the exper- 
imental phonon dispersion relations^ one can expect 
a large two-phonon density of states (DOS) for acous- 
tical modes spanning the range from about 80 to 160 
cm~^. So, decay of the zone-center optical phonons into 
two acoustical modes should be energetically possible. In 
view of the optical phonon softening and related strong 
anharmonicity (see also discussion below), an increase of 
the third-order coupling strength under pressure may be 
the dominant effect causing the line broadening. At this 
point we hesitate to elaborate further on this scenario, 
because we don't know the effect of pressure on the two- 
phonon DOS and, furthermore, pressure-induced changes 
in electronic structure should be considered as well. 

At ambient pressure, antimony is a semimetal with 
a negative indirect gap of '--^ 0.1 eV. The free carriers 
(density ^ 10^*'/cm'^) occupy states at Brillouin zone 
edges, forming pockets with small effective masses 
A coupling between the optical phonon modes and low- 
energy electronic excitations appears to be supported 
by a Raman study>i^ Now, the pseudogap of Sb is ex- 
pected to change significantly (towards filling up) under 
pressure, leading to a major change in the coupling of 
phonons to electronic states near the Fermi leveliS Fur- 
thermore, the response in the visible spectral range may 
be affected also, in particular the so-called E2 transition 
near 2.5 e V Is^ which at ambient pressure is in resonance 
with the exciting laser energy used in this study. So, we 
cannot neglect effects related to the penetration depth of 
the laser light. Furthermore, resonance conditions for 
Raman scattering may also change. 



III. PHONON FREQUENCY CALCULATIONS 

The phonon frequency calculations presented here 
are based on an evaluation of the total energy within 
the Density Functional Theory (DFT).^2iia We work in 
a plane-wave basis and use the Projector Augmented 
Waves (PAW) method^MS as implemented in the VASP 



TABLE IL Structural parameters of Sb-I obtained from DFT 
calculations using the LDA and GGA schemes. The values 
refer to the respective calculated zero-pressure volumes. The 
lattice parameters are for the non-primitive hexagonal cell, 
the volumes are for the primitive unit cell (2 atoms). The 
calculated values of bulk modulus Bo and its pressure deriva- 
tive B'q were obtained by fitting a third-order Birch equation 
of state to the calculated pressure- volume data (the attached 
errors reflect parameter correlations). Related experimental 
results are listed in the lower part of the table. 





a (A) 


c/a 


u 


Vo (A^) Bo (GPa) 


B'o 


LDA 


4.3003 


2.5406 
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" X-ray diffraction, Ref. [s^ 

X-ray diffraction, B'q fixed, Ref J3 

Ultrasonic wave velocities, Ref. |3^ adiabatic value 

Ultrasonic wave velocities, Ref. |38L adiabatic value 



codes4i The (pseudo-) potential^i*^ of Sb was available 
from G. Kresse and J. Joubert4^ Scalar relativistic cor- 
rections are taken into account indirectly, through the 
construction of the (pseudo-)potential. The expansion 
of wave-functions includes plane-waves with kinetic en- 
ergies up to the cutoff energy of — 215 eV. As 
antimony is metallic, we used for the Brillouin zone sam- 
pling a relatively dense uniform 14 x 14 x 14 mesh which 
corresponds to 280 k-points in the irreducible wedge4^ 
In order to deal with the partially occupied orbitals we 
employed the method of Methfessel and Paxton*^ (to 
the first order) with a smearing parameter tr = 0.1 eV. 
Whenever pressures are quoted in the context of calcu- 
lations, they were obtained from the DFT directly using 
the stress theorem. We have tested the LDA and GGA 
versionsiZ' of the (pseudo-)potential - essentially we re- 
peated all structural calculations (such as the E(y) and 
P{V) equations of state) in both LDA and GGA. Neither 
of the two approximations turned out to be convincingly 
superior to the other and most of the phonon calculations 
reported below relate to LDA unless stated otherwise. 

The phonon frequencies were calculated in the frozen 
phonon approach starting from optimized structural pa- 
rameters. Our optimized structural parameters at zero 
pressure (calculated) are given in Table ^ Results of 
structure optimizations at different volumes are discussed 
in the Appendix. 

The frozen-phonon approach consists in the evalua- 
tion of the total energy of the crystal with frozen- 
in atomic displacements. Small displacements v{l) = 
±(0,0, w) and v{2) — =f(0,0, u) in the Aig mode and 
v{l) = ±('i;,0,0) and v{2) = :^{v,0,0) in the Eg mode 
were applied to the two atoms of the basis, and the two 
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FIG. 6: Results of total energy calculations (LDA) for Sb in the Al structure. Left: Energy at different volumes as a function 
of phonon displacement for the Aig and Eg modes. In the order of decreasing volume, the calculated pressures are -1.68, -0.19, 
1.52, 4.14, and 5.46 GPa. In the Aig mode, the positive displacement corresponds to the shortening of the covalent Sb-Sb 
bond. Left: Restoring forces as a function of displacement for the Aig and Eg modes at two different volumes. In the case of 
the Eg mode, both the forces in x and z direction are given (Fy = at all volumes). For the Aig mode, there is no non-parallel 
force component. 



values obtained for Ai?*°* [viz. the Ai?*°* (outward) and 
Ai5'°' (inward)] were averaged. The eigenfrequency is 
then obtained from the expression for the energy of a 
harmonic oscillator 

t.E'°' = \mu?{\v{\)\^ + \v{n^\ , (1) 

where M is the atomic mass of Sb. As all the /S.E*-°^{v) 
variations suggested presence of anharmonic terms we 
chose the values of Vz/c = 0.004 (for Ai^) and Vx/a = 
0.012 (for Eg) at which the eventual quartic terms are 
negligible. The main anharmonic terms to take care of 
are the cubic ones, and they are eliminated by averaging 
the Ai;*°* (outward) and Ai;*°* (inward). 

Presence of the anharmonicity in the E(v) potential 
is clearly visible in Fig. El The scales of v are about 
the same for the Aig and Eg mode when measured in 
absolute displacements (0.1 A corresponds, at ambient 
pressure, to v^/c = 0.0089 or Vx/a = 0.023). We can 
see that, at low pressures, the energy variation in the Eg 
mode does not differ much from a parabolic shape: the 
main deviation shows as the slight asymmetry for positive 
and negative v. The degree of asymmetry remains about 
the same when pressure is increased. This behavior is in 
stark contrast with the one seen in the Aig mode where, 
at P = 0, the E{v) is essentially symmetrical at small 
displacements {v/c < 0.005) but its asymmetry rapidly 



increases under compression, to reach at = 53.14 A'^ 
(P « 5.5 GPa) an almost flat shape in one of the dis- 
placement directions. 

A closer inspection of the Aig potential energy curves 
reveals, at low pressures (P ~ 1.58 GPa), a small energy- 
barrier around Vz/c « —0.01; when pressure increases 
the barrier shifts to smaller values of ji'zl/c and becomes 
lower. An expansion of the unit cell, on the contrary, 
makes the barrier higher and shifts its position to larger 
values of the displacement (e.g. Vz — —0.015 c = 0.17 A 
at P = —1.68 GPa), where it is beyond the reach 
of the mean-square displacement at room temperature 
[vz ~ 0.09 A). It should be noted that the maximum 
of the barrier corresponds to the structural parameter 
u = 0.25. If we had plotted the results by measuring the 
displacement Vz relative to w = 0.25, we would have ob- 
tained symmetric potential curves similar to those shown 
in Refs.l(Sb) and|4i(Bi). 

More comprehensive information about the anhar- 
monicity is provided by the variation of forces with dis- 
placements, which is also plotted in Fig. El Now the 
harmonic regime is easily recognized by the linear vari- 
ation of the restoring force F against v. At all volumes 
close to ambient pressure the Aig mode is harmonic, until 
approximately the displacements of Vz/c = 0.004, which 
is the value we chose for the frozen phonon calculations. 
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to become conspicuously anharmonic a,t V = 54.18 
(P « 4.1 GPa). For the Eg mode, on the other hand, 
the Fx seems to deviate httle from the hnear (i.e. har- 
monic) behavior for v^/a < 0.02, even when pressure is 
increased. However, and perhaps to some surprise at the 
first sight, a component of force Fz appears in this 
mode. 

This somewhat unexpected emergence of the non- 
paraUel component Fz reminds us that a multitude of 
different anharmonic terms can appear in the expansion 
of E{v), and the presence of anharmonicity then can 
show up in two ways: (1) forces are not linear in v and 
(2) forces are not necessarily parallel to the displacement 
(see Ref. lifll Sect. III). It is the "mixed" VxVz terms in 

(vz) = Avl + Bvl + Cvlvx + Dvzvl + ... 
or (2) 
{vx) = Avl + Bvl + Cvlvz + Dvxvl + ... 

that are responsible for the "non-parallelism" of F and v. 
Figure IHl shows that, while the "parallel" term B seems 
to be significant in Aig, it is the large C which gives rise 
to the non-zero component Fz in Eg. Dealing only with 
the E{v) variations (Fig. |SJ) would, obviously, make us 
miss any mixed term in Eqs. 0. 

In the harmonic regime, using the Hellmann-Feynman 
forces merely offers a convenient, alternative way to ex- 
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FIG. 7: (Color online) Results of frozen-phonon calculations 
(LDA and GGA) for Sb in the Al structure: harmonic phonon 
frequencies at different unit cell volumes. The experimental 
data (dashed lines) are indicated for comparison. Experimen- 
tal pressure is converted to volume using a Birch equation of 
state with experimental parameters given in Table IITI 



plore the calculated Eiv) data. However, in an anhar- 
monic oscillator the knowledge of the variation F(v) pro- 
vides us with additional information on the mixed terms 
in Eqs. @. 

In the light of the above discussion, it does not come 
as a surprise finding out that the harmonic frequencies 
calculated from Eq. (1) agree rather imperfectly with 
the measured data shown in Fig.EI However, the overall 
trend of the variation is correctly reproduced. In par- 
ticular, the initial rate of the variation with volume is 
not far from the observed values and the larger nonlin- 
earity in the softening of Aig mode is reproduced also. 
The discrepancies remind us that the (calculated) har- 
monic frequencies are not the quantity that is measured 
in the spectrum of an anharmonic oscillator. Another 
contribution to the discrepancy is caused by a deviation 
of the optimized c/a ratios from experimental values (see 
Appendix) . 

An attempt to evaluate the frequency shifts or the line- 
broadenings in terms of some of the theories proposed in 
the pastS would lead us beyond the scope of this pa- 
per. In particular, one would have to assume that an- 
harmonicity is the only mechanism responsible for the w 
20 cm~^ shifts appearing in Fig.Q which is by no means 
certaini^ 



IV. CONCLUSIONS 

We have performed a high pressure Raman scattering 
study of Sb up to 16.6 GPa at room temperature. On 
the approach to the Sb-I to Sb-II phase transition near 
8 GPa, both the Aig and Eg zone-center phonon modes 
of the Al phase exhibit a pronounced softening, of the 
order -20%. Frozen-phonon calculations yield harmonic 
phonon frequencies in reasonable agreement with the ob- 
served mode softening. The two Raman lines of the A7 
phase show an unusually large broadening under pres- 
sure. The fact that the lifetime of phonons decreases 
with increasing pressure could become important if, sim- 
ilar to BifiSiSS the effect of pressure on coherent phonon 
generation in Sb were to become of interest. We report 
Raman-active modes of the phase Sb-II which has a com- 
plex guest-host structure. For the Sb-II phase, all the ob- 
served phonon features exhibit a hardening with increas- 
ing pressure, at least up to 16 GPa; a correlation with 
the negative pressure coefficient of the superconducting 
transition temperature of Sb-II-^ is apparent. 

Perhaps the most interesting question raised by the 
present study concerns the coupling of the zone-center 
optical phonons of the A7 phase of Sb to other phonons 
(e.g. anharmonic decay) and possibly electronic pro- 
cesses (electron-phonon coupling). It is the interplay 
between mode softening, electronic structure changes at 
the Fermi level, and shifts of optical interband tran- 
sitions (i.e. resonances and penetration depth) which 
leads to a more complex scenario of phonon decay pro- 
cesses compared to semiconducting materials. Our low- 



8 



temperature measurements of the Raman line widths in- 
dicate that phonon-related processes are a major cause of 
the pressure-induced broadening. So, a first step towards 
gaining more insight would be to study the effect of pres- 
sure on the phonon dispersion relations, specifically the 
acoustic branches. 



APPENDIX: A7 STRUCTURE OPTIMIZATION 

The literature about different phases of anti- 
mony includes several studies based on first-principles 
calculations. An interesting problem had been the 
question of the existence of a cubic primitive (cPl) struc- 
ture of Sb at pressures between 6 and 9 GPa. The cPl 
structure can be viewed as a special case of A7 with 
u — 0.25 and c/a = ^/6. First suggested experimentally 
in Ref. Is^ the reality of the cPl phase was no t co n- 
firmed in later diffraction studies of Sb, e.g. Refs. I2l6l l3. 
Since then the A7 — > Sb II/IV transition is believed to be 
the first in the sequence of pressure-induced transitions 
of Sb, and the 'cubic-primitive' subject has occasionally 
appeared only in theoretical studies which slowly "con- 
verged" to elimination, on energetic grounds, of the cPl 
phase from the list of stable phases of antimony. 

Our calculated results for the optimized structural 
properties of Sb-A7 as a function of unit cell volume, 
in part sampled at very small volume intervals of 0.1 A'^, 
are summarized in Fig. |H1 Let's first consider the LDA 
results. Starting from large volume, the calculated vari- 
ation P{V) [Fig. IHIa)] is smooth first. When passmg 

V 52.9 A3 a kink appears in P{V), followed by narrow 
interval of smooth variation, until V = 52.1 A^. The 
kink in P{V) relates to the evolution of the structural 
parameters c/a and u [Figs. |Hl[b,c)]. At the approach to 

V = 52.9 A^ the smooth decrease of c/a changes slope 
and c/a takes, practically discontinuously, values situ- 
ated close to the cubic value c/a — V6. At the same vol- 
ume, a nearly abrupt variation towards the cubic limit 
u — 0.25 is also observed in the behavior of u. How- 
ever, the cubic limit of u is only reached by a second 
'jump' near V ^ 52.1 A^. At all volumes V < 52.1 A^, 
the energy-minimization algorithms, set to conserve the 
i?3m symmetry, steer the iteration obstinately towards 
exactly u = 0.25, but with c/a significantly larger again 
than the cubic value c/a = V6. When increasing the 
volume from below to above V = 52.1 A^, the optimiza- 
tion procedure tends to stay trapped at u — 0.25 until at 
about y = 54 A'^ it falls back to the curve obtained for 
decreasing volume. 

The corresponding E{V) calculations (not shown here) 
tell us that the range V G (52.1 ; 52.9) A^ (pLDA between 
5.7 and 6.6 GPa) should be considered as the region of 
co-existence of two A7 variants, because the total energy 
differences are smaller than 2 meV/cell and the specific 
volumes are almost identical. The third A7 variant is 
energetically favored for F < 52.1 A^ (pLDA > g g Qp^y 
However, at these volumes the A7 — > Sb-II transition (not 
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FIG. 8: (Color online) Calculated structural properties of 
Sb in the A7 structure as a function of rhombohedral unit 
cell volume (2 atoms per cell): (a) equation of state P{V); 
(b) axial ratio c/a; (c) positional parameter u. In the LDA 
approximation, a sequence of two structural transitions to 
variants of the A7 structure (see text) is predicted between 
V = 52.9 and 52.1 A"^ (marked by dashed vertical line). In 
GGA, the corresponding volume range is larger. Closed cir- 
cles and open squares represent the tabulated experimental 
resultsi-*- The dashed vertical line at 52.1 A'^ coincides with 
the smallest volume for which the pure A7 phase has been 
observed experimentally at 300 



investigated here) is even more favored; the experiments^ 
situate the appearance of the latter phase at V — 52.2 ± 
0.02 A'^. Hence, the possibility of observing that third 
A7 variant can be ruled out. 

The peculiar variations of c/a and u at V — 52.9 A^ 
indicate a transition to a pseudo-cubic A7 variant where 
only the metric approaches the cubic value (c/a = v^)- 
The symmetry of the structure remains R3m. The sec- 
ond discontinuity appearing in u and c/a at ^ = 52.1 A'^ 
suggests a transition towards a a third A7 variant where 
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the behavior of cja and u is reversed; the 'cubic' value of 
u = 0.25 exactly, constant and independent of volume, 
and cja > V6 and varying. So, allowing for the two 
structural degrees of freedom of Rim, our calculations 
never converged towards a stable cPl structure of Sb. 

The non-smoothness of the c/a{V) and u{V) varia- 
tions, occurring at essentially the same two volumes as 
marked by vertical lines in Fig. |H1 was spotted in the 
context of comparisons of LDA vs. GGA calculations re- 
ported in Ref. jsTf The details of the variations were not 
pursued throughout the "instability region" , and Ref. 
concluded at a phase transition to the cPl structure at 
V <b2 K? (incorrect) , but found out that "the axial ratio 
converges to c/ a^/& before the parameter u converges to 
u = 1/4" (correct). Another, more recent, calculations- 
concluded that the A7 Sb-II transition takes place at 
P = 8.6 GPa, i.e. before the (hypothetical) A7 cPl 
transformation might occur at the (calculated) pressure 
of P = 10.7 GPa. Apparently, non-smooth variations of 
the A7 structure were not examined in detail in Ref. 12^ . 

So as to estimate the uncertainties of our above proce- 
dures, we briefly compare the LDA calculations with the 
ones performed using the GGA. Within that approxima- 
tion a similar behavior was obtained but the "coexistence 
region" of A7 variants is larger, i.e. V G (52.2 ; 53.8) A'^. 
A sensible estimate is thus to situate the region of in- 
stabilities, which is to be compared to experiment, as 
the average of the two results: V e (52.15; 53.35) A^. 
Using then also for the V ^ P conversion the aver- 
age of the LDA- and GGA-limits for P{V), we obtain 
a lower "critical" pressure for the instability region as 
P « 7.3 GPa. This is lower than the experimentally ob- 
served appearance of the Sb-II structure (8 to 8.5 GPa 
on upstroke). But, taking into account that the Sb-I to 
Sb-II transition exhibits a hysteresis of about 1 GPa at 
room temperaturef2i2& the calculated critical pressure is 
only slightly less than the equilibrium transition pressure 
(roughly estimated as being the average pressure value 
for the upstroke and downstroke). So, there is not much 
room for actually finding an instability of the A7 phase 
preceding the Sb-I to Sb-II/IV transition. 

In reports on early x-ray diffraction studies performed 
under non-hydrostatic conditions, the authors uphold 
to have observed the sluggish transformation to an in- 
termediate primitive cubic structure of antimony in the 
pressure interval P e (7.0; 8.5 — 9.0) GPa or at least pro- 
nounced anomalies in the evolution of the A7 lattice pa- 
rameters with pressure. The findings of Ref. Is^ derived 
from low-resolution diffraction diagrams, have not been 



confirmed in subsequent diffraction studies. Also, based 
on discontinuities in the pressure variation of the elec- 
trical resistivity of Sb seen at P = 7.8 GPa (upstroke) 
and 6.7 GPa (downstroke) the existence of only one 
structural transition (A7 to Sb-II) is proposed in the in- 
terval considered here. Nevertheless, the parallel between 
our calculated results and the experimental observations 
of Ref. .52 is. perhaps, not accidental. It is conceivable 
that, in the presence of a substantial uniaxial stress com- 
ponent along the C3 axis, a precursor transition can be 
induced. In this context it should be noted that, at fixed 
volume, the calculations yield optimized c/a ratios which 
are smaller than the experimental values [cf . Fig. IH^b)] . 
So, the condition of zero deviatoric stresses in the cal- 
culations corresponds to a strained configuration when 
compared to experimental data. With attention now be- 
ing drawn to the importance of uniaxial stress, we arrive 
at what could be expected from intuition: the Peierls-like 
distortion of Sb-I, being a uniaxial effect, may depend 
sensitively on uniaxial stress along C3, with associated 
structural variations that are not observed under truly 
hydrostatic conditions. 

To summarize, total energy calculations constrained 
to space group R3m point to the possible existence 
of an A7 to A7 transition near the borderline where 
antimony transforms to the higher-pressure Sb-II/IV 
phases. We do not rule out other modes of instability 
of the A7 structure that require a space group lower in 
symmetry (monoclinic) than R3m. So, a more general 
conclusion would be that the 'normal' A7 structure 
of Sb may exhibit a structural instability before it 
transforms to the Sb-II/IV phases. Some tuning of the 
phase behavior through uniaxial stress may be possible. 
When comparing the calculations with experiment, an 
inherent uncertainty concerns the role played by a finite 
temperature. The present as well as other ab initio 
calculations of Sb were performed in the static lattice 
limit. 
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